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Solid-fluid transition in a granular shear flow 
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The rheology of a granular shear flow is studied in a quasi- 2d rotating cylinder. Measurements 
are carried out near the midpoint along the length of the surface flowing layer where the flow is 
steady and non-accelerating. Streakline photography and image analysis are used to obtain particle 
velocities and positions. Different particle sizes and rotational speeds are considered. We find a 
sharp transition in the apparent viscosity (77) variation with rms velocity (w). In the fluid- like 
region above the depth corresponding to the transition point (higher rms velocities) there is a rapid 
increase in viscosity with decreasing rms velocity. Below the transition depth we flnd 77 oc u~^'^ for 
aU the different cases studied and the material approaches an amorphous solid-like state deep in 
the layer. The velocity distribution is Maxwellian above the transition point and a Poisson velocity 
distribution is obtained deep in the layer. The observed transition appears to be analogous to a 
glass transition. 

PACS numbers: 45.70.-n, 45.70.Mg, 83.80.Fg 



Granular materials are known to exist in solid-like and 
fluid- like states . Physical understanding of the flow of 
granular materials has thus developed along two major 
themes based on the flow regime 3 • In rapid flow — 
fluid-like — regime, both theory and experimental analy- 
sis are generally cast in the framework of the kinetic the- 
ory In contrast, the slow flow — solid-like — regime 
is most commonly described using the tools of soil me- 
chanics andplasticity theory 4] and recently by analogy 
to glasses S • ^ '^'^y difference between solid-like 
and fluid-like states, which appears to be emerging from 
recent studies, is the definition of a temperature in the 
two states. In the fluid-like state the granular temper- 
ature is defined as the kinetic energy of velocity fluctu- 
ations in analogy with the kinetic theory of gases. In 
contrast, for dense granular flows (solid-like state), nu- 
merical 111 and experimental studies show the validity 
of the fluctuation-dissipation theorem (the mean-square 
displacement of a particle is proportional to time for a 
constant applied force). This allows for the definition of 
a new temperature 7]. These two approaches have no 
well-understood region of overlap. Given the qualitative 
differences between the fluid and solid states, a question 
that has been open for some time relates to the criterion 
for transition between fluid and solid states of granular 
materials. We focus on this question. 

There are few studies which focus on the transition 
between the solid and fluid states. Metcalfe et al. 
studied solid to fluid transition in a horizontally vibrated 
container of beads. They observed hysteresis in the tran- 
sition which was well-predicted by a dry friction model 
in which the friction coefficient varies smoothly between 
a dynamic and static value. A fluid-solid transition was 
also observed by D'Anna and Gremaud ,5(| for vertically 
vibrated particles. They found that the vibrating granu- 
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FIG. 1: (a) Schematic of the rotating cylinder geometry show- 
ing the flowing layer and the co-ordinate system employed, 
(b) A typical image showing the portion of the flowing layer 
at the center of the cylinder, with streaks of various lengths 
formed across the layer for a shutter speed of 1/250 s. (c) 
Magnified image of the rectangular region marked in (b) show- 
ing streaks generated by three different particles, (d) The 
same three streaks with the fitted intensity function. The 
optimization yields length (2Z), width (2u;), orientation angle 
(9) and position {x, y) for every streak. 



lar medium first undergoes a sharp transition to a super- 
cooled liquid and then gradually achieves a solid-like 
state on reducing the intensity of vibration. This pro- 
cess follows a modified Vogel-Fulcher-Tammann (VFT) 
model which is typical of fragile glasses. 

Coexisting solid and fluid phases have been studied 
primarily in the context of surface flows; these comprise 
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FIG. 2: Mean velocity {{vx)), rms velocity (u), number density (n), shear rate (7) and shear stress {jyx) across the flowing 
layer for 2 mm steel balls rotated at 3 rpm. The solid line in (a) shows a linear fit to the mean velocity profile. The dashed 
line denotes the transition point in all the figures. 



a layer of fluid like flow on a fixed bed of the same ma- 
terial. Examples of surface flows which have been well- 
studied include heap flows [il l fia. |l ^ [la lia llTl 
and rotating cylinder flows psl Il9l l20l l2ll 123. 123 . l24l| . 

Remarkably simple theories describe the coexistence be- 
tween the solid and fluid. In the simplest versions, the 
local melting and freezing is determined by the local an- 
gle of the solid-fluid interface: If the local angle is greater 
than a "neutral" angle, the solid melts (the heap erodes) 
so as to reduce the angle and vice versa 0, . Con- 
tinuum models based on Coulombic friction models for 
the solid region and simple rheolqgical models for fluid 
region predict a similar behaviour |l4l IT^ . Experimental 
studies show that the models based on local angle based 
melting/freezing give good predictions T^. 

An assumption in the coarse-grained models described 
above is the existence of an interface between the solid 
and fluid. However, the recent work of Komatsu et al. 
[TEl l indicates that the surface flow on a heap decays 
smoothly with depth and generates motion deep within 
the heap. The interface between the solid and fluid re- 
gions is thus not well-defined when considered at a parti- 
cle length scale. Studies of the velocity profile in rotating 
cylinders flow also confirm this picture of a smooth de- 
cay of the velocity into the bed, rather than an abrupt 
change at an interface [2^[23l |. 

The objective of the present work is to gain an insight 
into the fluid-solid transition in granular flows focussing 
on a system where both fluid-like and solid-like regions 
coexist. We find a well-defined transition point in the 
system which demarcates two distinct flow regions. The 
behaviour in each region is characterized and results in- 
dicate that the transition has some characteristics of a 
glass transition. 

Experiments are carried out in quasi-2d aluminium 



cylinders (length 1 or 2 cm) of radius 16 cm (Fig. ^). 
The end walls are made of glass and a computer con- 
trolled stepper motor with a sufficiently small step is used 
to rotate the cylinders. Monodisperse, spherical, shiny 
stainless steel balls of three different sizes with diameters 
(d) 1, 2, and 3 mm are used in the experiments. Cylin- 
der rotation in the rolling flow regime (rotational speed, 
Lo — 2 — % rpm) produces a thin flowing surface layer 
and measurements are made near the center of the cylin- 
der where the layer thickness is maximum and the flow 
is non-accelerating. The particles are heavy enough and 
conductive so that charge effects are negligible. The ex- 
periments are carried out with 50% of the cylinder filled 
with particles. 

The motion of the particles is captured by taking high 
resolution images using a digital camera (Nikon Coolpix 
5000) in the presence of an incident beam of light. The 
size of the recorded region is 2560 x 1920 pixels, with 
one pixel corresponding to 0.016 — 0.03 mm depending 
on the distance of the camera to the cylinder. The point 
source of light is directed near parallel to the end wall 
of the cylinder so as to illuminate only the front layer 
of the flowing particles. Each moving particle generates 
a streak of definite length depending on its speed and 
the shutter speed of the camera. Images are taken for a 
range of camera shutter speeds (1/15 — 1/2000 s) so as to 
account for the varying velocity across the flowing layer. 
This gives streaks, which are adequately long for analysis 
but not so long as to overlap with other streaks, in each 
the part of the flowing layer (Fig. Qd). Two hundred 
images are taken for each shutter speed with an overall 
of two thousand images combined over different shutter 
speeds. Due to the time interval between photographs, 
each experiment typically takes 500 cylinder revolutions. 

A parameterized intensity function corresponding to a 
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stretched Gaussian function is then fitted to the intensity 
values of the streak pixels (and an immediate neighbour- 
hood) in the image by means of a standard non-linear op- 
timization technique. The fitting yields the length (2^), 
width (2w), orientation angle (9) and the position {x,y) 
for the streak (Fig.^). The analysis technique was cal- 
ibrated by carrying out experiments for a single particle 
glued to the cylinder end plate. The error in velocity 
measurements was found to be less than than 3%. For 
analysis, the flowing layer region is divided into bins of 
width equal to the particle diameter and length 20 mm 
parallel to the flowing layer. The components of the mean 
and the root mean square (rms) velocities for each bin are 
calculated by averaging over all streaks in a bin. 

Fig. El gives the variation of system variables with 
depth (y) in the flowing layer for 2 mm particles. The 
mean velocity {{vx)) profile is smooth and shows 3 re- 
gions; a near-linear middle region, an exponentially de- 
caying region at the bottom and a flattened region near 
the top. Similar profiles have been reported in several 
previous studies [Hlilillilllllillillll. The flat- 
tened upper region which is the very low density region 
corresponding to saltating particles is not seen in some 
studies. The magnitudes of the rms velocities in the two 
directions are different: the component in the flow di- 
rection (x-component) is about 10% higher than that in 
direction perpendicular to the flow (y-component). The 
profile of the total magnitude of the fluctuating veloc- 
ities (u) is shown in Fig. ^p. The rms velocity pro- 
file shows two distinct regions: a relatively slow varia- 
tion near the free surface followed by a sharper decrease 
deeper in the bed. Fig. |3i shows the shear rate varia- 
tion in the layer, obtained by numerical differentiation of 
the data in Fig. The shear rate increases to a maxi- 
mum value (corresponding to the inflection point in the 
velocity field) and then decreases. We note that an oscil- 
lating shear rate profile is obtained if a smaller bin size is 
used as reported previously 0i US ■ The transition 
point in the fluctuation velocity profile coincides with the 
maximum in the shear rate. The areal number density 
is almost constant throughout the fiowing layer with a 
rapid decrease near the free surface as shown in Fig. O;. 

For non-accelerating fiows a force balance yields the 
shear stress as dTxy/dy — pgsm/S « pb{n/nb)gsinP 
where g is acceleration due to gravity, P is the angle of 
repose, n is the number density in the fiowing layer, nt, is 
the number density in the rotating packed bed and pb is 
the bulk density of the rotating packed bed. We neglect 
the contribution of wall friction in the estimate. This 
could become significant deeper in the bed 28] . However, 
based on the method of Taberlet et al. we find that 
the contribution is about 10% of the total stress in the 
cases studied and does not qualitatively affect the results. 
Upon integration we obtain r^y — {pb/nb)g sin fi ndy 
assuming Txy = at the free surface. The shear stress 
shows a near linear increase with depth (Fig. |3i). The 
results of Fig.|21indicate that there is a qualitative change 
in the rheology at the transition point (j/c)- Above yc, the 
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FIG. 3: Apparent viscosity (77) variation with rms velocity 
(it) for all the cases studied. Filled symbols represent the 
points at and above the transition velocity (uc). (a) The 
solid line at the top represents a linear fit (slope ~ 1.5) to the 
data below the transition point for all the cases studied, (b) 
The data is scaled using particle mass (m), particle diameter 
(d) and average shear rate obtained by a linear fit to each 
corresponding velocity profile. 



shear stress increases with shear rate which is typical of 
fluids. However, below yc the shear stress increases while 
the shear rate decreases. This implies that the viscosity 
increases sharply with depth below yc, even though the 
number density is nearly constant. 

Fig. shows the variation of the apparent viscosity 
[f] ~ Txy/'y) with the rms velocity (w), considering only 
the region of constant density. The data points at and 
above the transition point (y > yc) are plotted as filled 
symbols. There is a sharp transition in this case as well 
and the transition occurs at the same value of u as in 
Fig. [2b- In the region near the free surface there is a 
rapid increase in viscosity with decreasing rms velocity 
whereas in the region approaching the fixed bed there 
is a much slower power law increase with an exponent 
n —1.5 for all the cases studied. The data scales with 
javg and d to fall on a single curve (Fig. ^p), and the 
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FIG. 4: Distributions of velocity in y direction at various lo- 
cations across the layer for 2 mm steel balls rotated at 3 rpm. 
The dashed line with filled circles represents the transition 
point, the solid lines represent the region above the transi- 
tion point and all the dotted lines represent the region below 
the transition point. The thick black line represents a fitted 
parabola to the curves at and above the transition point. 



transition occurs at Uc ~ \.2'^avgd for all cases. Here 
iavg is the shear rate obtained by fitting a straight line 
to the linear portion of the velocity profile (Fig.|3i.). 



Fig. 0] shows the distributions of the y direction veloc- 
ities at different locations in the layer. The Vy distribu- 
tion is Gaussian for all points above the transition point 
{y ^ Dc) and it gradually evolves to a Poisson distribution 
as we go deeper into the bed. A Poisson velocity distri- 
bution was also found by Mueth |23| for Couette flow 
of a dense granular material. The Vx distribution (not 
shown) is Gaussian above the transition point. However, 
below the transition point the behaviour is complex and 
bimodal distributions are obtained. 

The results presented show a sharp transition between 
two flow regimes. In the upper region near the free sur- 
face, the behaviour is fluid-like and the velocity distribu- 
tions are Maxwellian. Below the transition point the ma- 
terial appears to be an amorphous soft solid, increasing 
in strength with depth in the layer. The transition to this 
solid-like regime occurs at a relatively large rms velocity 
(w 0.1 m/s). We conjecture that the sharp transition oc- 
curs because of the formation of a percolated network of 
particles in extended contact with each other. This is in 
contrast to the fluid-like regime where the particles inter- 
act through collisions. The contact network coexists with 
fluid-like domains and the fraction of particles which are 
part of the network increase with depth. The transition 
appears to be analogous to a glass transition. New rhe- 
ological models may be required for granular flows that 
span these two regimes. 
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